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AFOSR.TH.  8  7-116? 


Pinal  Report 

Grant  No.  AFOSR-83-348 

Studies  of  the  structural  Dynamic  Behavior  of  Satellite  Antenna  Systems 


I  Introduction 

The  "hoop -column"  or  "hoop -maypole"  type  of  antenna  (Pig.  1.)  is  typical  of 
many  of  the  space  structures  being  designed  and  anticipated  for  use  in  near-earth 
orbit.  The  natural  modes  and  frequencies  of  such  sparse  structures,  which  are 
exceptionally  large  and  light  (i.e. ,  flimsy  by  earthly  standards),  must  be  known  with 
considerable  accuracy  to  insure  that  corrective  impulses  (or  other  control  forces  and 
moments),  as  applied  by  vernier  thrusters  or  control  gyros,  will  not  inadvertently 
trigger  unstable  closed -loop  oscillations. 

Motivation  for  the  subject  studies  were  two  foldi  first,  to  gain  some  insight  as  to 
the  changes  in  natural  modes  and  frequencies  with  large  variations  in  design 

parameters;  and  second,  to  develop  a  transfer  matrix  method  of  analysis  as  an 

alternative  to  finite  element  methods  (PPM)  which  are  presently  the  primary  analysis 
tool  in  dealing  with  the  dynamic  response  of  complex  space  structures . 
Straightforward  applications  of  powerful  analysis  approaches  such  as  the  finite  element 
method  leave  much  to  be  desired .  The  results  are  usually  of  such  volume  and 

complexity  as  to  obscure  the  nature  of  the  responses  and  hinder  the  proper  choice  of 

design  changes.  A  structural  breakdown  sufficient  for  accuracy  in  such  analyses  is 
also  costly  in  manhours  and  computer  time.  The  transfer  matrix  approach  promises 
efficiency  at  a  level  which  is  appropriate  for  a  practical  design  tool;  at  the  same 
time  its  modular  nature  may  increase  insights  and  understanding  regarding  the  dynamic 
behavior  of  flexible  space  structures. 

The  work  reported  here  was  conducted  in  several  phases;  (a)  formulation  of  the 
complete  analysis,  (b)  two-dimensional  trend  studies  of  flexible  beam-columns  from 
which  a  rigid  bar  is  suspended  by  cables,  (c)  two  and  three  -dimensional  trend  studies 
of  "T"  and  "H"  structures  typical  of  various  combinations  of  substructure  consisting  of 
mast,  feed  assemblies  and  solar  panels,  (d)  trend  studies  of  the  modes  and 
frequencies  of  planar,  polygonal,  cable -stiffened  "hoop"  structures  and  (e) 
investigations  into  the  numerical  difficulties  enountered.  Pig.  2  (from  Ref.  l)  is 
representative  of  the  kinds  of  variation  likely  to  be  encountered  in  design  studies. 

Cl  Results 

Phase  A  (Theoretical  Development) 

The  first  phase  resulted  in  a  transfer  matrix  formulation  in  which  bending  in  two  — - 
normal  directions,  compression/extension,  and  torsion  are  all  represented  in  the  major  V 
structural  components  of  the  system;  namely  feed  assembly,  central  column  (or  mast),  □ 
hoop  and  solar  panels.  Cables  are  assumed  to  run  from  the  vertices  formed  by  the  n 
ends  of  each  hoop  element  to  three  different  longitudinal  stations  on  the  central  mast; 
they  are  treated  as  massless  springs  in  which  there  are  (large)  preloads.  The 
portion  of  the  mast  between  hoop-cable  attachments  carries  a  (large)  steady,  ... 
compression  load,  as  do  the  hoop's  polygonal  sides.  The  feed  assemblies  are  treated 
as  attached  to  the  mast  at  an  angle  of  90°,  but  also  have  an  arbitrarily  large  "bend”  - 
built  into  them,  to  account  for  "aiming"  their  E-M  radiation  properly  at  the  reflecting  °dos 
surface.  The  solar  arrays  are  treated  similarly,  but  here  an  arbitrattly  large,  7 


discreet  "twist"  is  introduced  to  account  for  their  being  oriented  toward  the  sun. 
Shear  deflections  are  accounted  for,  as  well  as  mass  -offsets ,  which  have  the 
consequence  of  eliminating  symmetry. 

The  full  analysis,  then,  involves  12  x  12  transfer  matrices,  which  successively 
premultiplied,  as  unknown  state  vectors  are  eliminated ,  lead  to  a  6  x  6  determinant 
whoee  value  should  be  zero  for  the  correct  choice  of  a  trial  frequency.  The  solution 
procedure,  therefore,  is  to  -  by  trial  and  error  -  find  the  frequencies  which  make 
this  determinant  zero  (Ref.  2). 

All  the  results  reported  below  for  the  phases  of  the  project  subsequent  to  this 
one  used  special,  more  limited  cases  of  the  theory  developed  in  Phase  A. 

a  (Two  Dimensional  Dynamics  of  a  Rigid  Bar  cabled  to  a  Flexible  Mast) 

For  this  phase  of  the  research,  more  hilly  reported  in  Ref.  3,  the  baseline 
antenna  configuration  was  assumed  to  have  the  characteristics  given  in  Table  1,  and 
no  steady  load  in  the  cables.  The  mathematical  model  for  the  flexible  mast  had  five 
discreet  masses  outside  the  cable  attachment  points  and  9  inside  them  for  a  total  of 
19  masses,  in  all  mast  length  variations,  the  mass  per  unit  length  was  held  constant; 
in  hoop  (rigicl  bar)  length  variations,  the  mass  per  unit  length  was  held  constant  but 
its  mass  moment  of  inertia  was  assumed  to  vary  as  that  of  a  rigid  circular  hoop. 
Changes  in  natural  frequencies  as  the  mast  length  varies  are  shown  in  Figs.  3a  and  b. 
These  configurations  are  symmetric  above  and  below  and  hold  the  portion  of  mast 
length  between  cable  attachments  constant.  The  odd  numbered  bending  modes  are 
symmetric  modes;  those  even  numbered  are  antisymmetric  about  the  configuration's 
center.  The  half -structure  analysis,  using  appropriate  boundary  conditions  for 
symmetry  and  antisymmetry  gave  results  identical  (as  expected)  to  those  of  the  full 
structure  and  resulted  in  30X  computer  running  time  savings.  Proper  boundary 
conditions  for  the  half  structure  are; 

For  the  symmetric  case: 
at  the  mast  center  - 

lateral  bending  slope  =  0 
lateral  shear  force  -  0 
axial  displacement  =•  o 

for  the  half  bar  - 

rotation  -  0 

x  and  y  components  of  cable  force  equal  the  D'Alembert 
forces  due  to  those  components  of  the  (one  half)  bar's 
translation 

For  the  antisymmetric  case: 
at  the  mast  center  - 

lateral  displacement  =  o 
bending  moment  -  0 

axial  force  =  0 

for  the  half  bar  - 

bar  center  displacement  3  0 

the  cable  force  components  perpendicular  to  the  bar  times  the 
bar  length  equal  the  O' Alembert  moment  due  to  the  (one  half) 
bar  rotation. 


3. 


When  the  length  of  the  mast  above  the  cable  attachment  points  is  varied  (ie  the 
length  of  the  "feed  mast"),  keeping  its  mass /length  and  all  other  mast  properties 
corartant,  one  obtains  the  trends  shown  in  Pigs.  4a  and  5.  This  is  one  of  three  kinds 
of  asymmetry  ( lengthwise )  examined  in  Ref.  3.  Among  other  trends  similarly 

calculated  were  the  effect  of  varying  bending  stiffness  outside  the  cable  attachment 
points  relative  to  that  inside,  the  rigid  hoop  (bar)  length  and  cable  attachment 

angles.  Trends  where  the  large  steady  axial  load  in  the  mast  between  cable 
attachments  was  added  and  varied  relative  to  its  Euler  buckling  load  were  also 

obtained.  While  frequencies  dropped,  as  expected,  no  significant  mode  shape  changes 
occurred  to  within  a  factor  of  2  of  the  buckling  load.  To  assess  the  effects  of 

moving  a  suspended  hoop  (bar)  from  its  centrally  attached  position,  all  the  mast 

properties  were  held  constant,  and  the  cable  attachment  points  moved  axially  to 

various  offsets  from  the  initially  symmetric  position.  Cable  angles,  bar(hoop)  length, 
etc  all  remained  unchanged .  The  resulting  variations  are  shown  in  Pig.  5a  and  b. 
(in  these  figures  Pa  =  axial  compressive  load  in  the  length  of  Central  column  between 

cable  attachments ,  and  Pg  =  the  critical  buckling  load  in  that  member.)  The 

existence  of  steady  axial  loads  induced  by  cable  tension  is  more  influential  as 

symmetry  is  lost  for  some  modes,  and  less  influential  in  others.  A  total  of  32 
different  configurations  are  considered  in  Ref.  3  with  up  to  10  natural  modes  and 
frequencies  given  for  each  configuration. 

This  simple  model  provides  a  useful  reminder  that  for  such  systems,  the  classical 
count  of  nodes  (points  of  zero  deflection)  is  only  a  reliable  indicator  of  mortal 
number  if  one  can  make  a  separation  between  modes  primarily  involving  the  continuous 
structure  and  those  involving  soft -sprung  substructures.  This  is  clearly  shown  in  Pigs. 
6a,  b  and  c. 

Phase  C  (Two  and  Three  Dimensional  Dynamics  of  "T"  and  "H"  substructures) 

This  phase  continued  studies  of  component  substructure  dynamic  behavior t  first, 
the  "T"  formed  by  the  feed  assemblies  or  the  molar  panels  mounted  on  a  central  mast, 
and  then  the  "H"  (really  on  its  side)  formed  by  the  feed  assemblies,  mast  and  solar 
arrays.  These  were  first  exaadned  as  two-dimensional  structures.  Tables  2,  3  and  4 
list  the  properties  used,  which  were  arrived  at  on  the  basis  of  some  known  distensions 
and  materials  and  assumptions  intended  to  be  comps  table  with  those  of  the  hoop 
maypole  antenna  described  in  Ref.  4.  Both  "T"  and  "H"  substructures  can  be 
symmetric,  and  for  such  cases  the  half -"T~  and  half -"H"  were  analyzed  and  identical 
results  obtained  with  check  calculations  using  the  full  structure  (to  within  four  places 
of  Un  value).  Proper  boundary  conditions  for  the  "T~  and  the  "H"  where  the  full 
symmetric  structure  is  split,  ie  along  the  center  of  the  mast,  are: 

for  the  symmetric  case  - 

Transverse  Shear  Pores  3  o 

Bending  Moment  3  0 

Bending  Slope  3  0 

Lateral  Displacement  3  o 

for  the  antisymmetric  case  - 

Axial  Displacement  3  o 

Axial  Pores  3  o 


Boundary  condition*  were  zero  fore*  and  moment  at  the  fra*  end*  of  both  structures; 
the  baa*  of  the  “T",  however,  was  assumed  to  be  cantilevered ,  ie  aero  deflection* 
and  rotation*  at  that  point. 

Chance*  in  frequency  and  node  shape*  were  calculated  as  the  length*  of  feed 
assemblies  and  solar  panels  were  varied ,  and  for  these  cases  of  the  ~T~ 
substructures,  some  general  characteristics  can  be  cited. 

o  Solar  panels  are  so  flexible  in  transverse  bending  motion,  that  no  axial 
extensions  of  the  mast  are  perceptable  in  symmetric  transverse  bending  motions  up  to 
the  sixth  mod*. 

o  The  Solar  Panel  "Base's"  transverse  bending  rigidity  (or  lack  of  it)  is  such 
that  the  "continuity  of  slopes"  required  of  symmetric  mod*  shapes  appears  to  be 
violated  (and  is  not) . 

o  Solar  Panel  Mass/Stiffness  distributions  are  such  as  to  result  in  unusually 
low  vibratory  deflections  near  the  free  ends  in  the  first  six  modes. 

o  The  order  of  feed  assembly  /upper  mast  modes  (l**  a  2nd  antisymmetric ,  3rd 
symmetric,  4*1  a  5th  antisymmetric,  6"1  symmetric)  is  unchanged  as  feed  assembly 
length  is  varied  from  1/2  to  3  times  the  base  case  (with  constant  mass/unit  length) . 

°  The  6th  feed  aaseably/upper  mast  mode  (2nd  symmetric)  shows  significant 

axial  mast  motion {  the  smaller  the  feed  assembly,  the  greater  this  motion. 

o  Reduction  in  CPU  time  using  the  half  "T"  compared  to  the  full  T"  was 
26.51. 


variations  of  the  symmetric  "H"  substructure  included  solar  panel  length  increases 
(outward)  to  130*  of  its  base  caset  increasing  solar  panel  length  inward  ie  shortening 
attachsmnt  and  foundation  arms  comsmnm  irately )  to  130*  j  increasing  mast  length  to 
163*,  increasing  mast  stiffness  by  factors  up  to  8,  increasing  feed  assembly  stiffness 
by  factors  up  to  6  and  length  (hence  mass)  changes  of  from  1/2  to  3  timee  the  bass 
case.  Several  general  conclusions  can  be  drawn  as  regards  the  first  six  natural 
n»des  and  frequencies. 

o  they  appear  to  be  dominated  by  mast  and  feed  assembly  mams  and  by  solar 
panel  (plus  attachment)  flexibility. 

o  the  order  of  modes,  first  antisymmetric,  then  alternating  to  the  sixth  (which 
is  symmetric)  is  unchanged  throughout  the  complete  rang*  of  parameter  variations. 

o  consistent  with  the  alternating  (antisymmetric/symmetric)  nature  of  the  first 
six  modes  and  the  mass  dominance  of  the  mast  and  feed  assembly,  natural  frequencies 
are  closely  "paired".  That  is,  each  pair  of  symmetric  and  antisymmetric  modes  are 
doss  in  frequency,  with  the  usual  sort  of  frequency  separation  (factors  of  2  to  3) 
between  one  pair  and  another. 

o  the  effect  of  the  mast  stiffness  changes  were  not  reflected  in  natural 

frequencies  (to  four  places)  and  were  imperceptabl*  in  the  node  shapes. 

o  the  effect  of  feed  assembly  stiffness  and  length  (and  therefore  mass 

changes  were  perceptable,  but  very  slight. 


> 


Additional  trends  war*  obtained  for  changes  which  destroyed  the  symmetry  of  the  "T" 
and  "H"  aubatructurea .  For  the  "T*  casee  examined,  auch  changes  caused  the  two 
aides  (feed  assemblies  or  solar  arrays)  to  behave  much  more  independently.  For  the 
“B"  cases,  the  symmetry /antisymmetry  mode -type  alternation  was  broken  and  so  was  the 
frequency  "pairing". 

Complete  results  are  given  in  Ref.  5,  where  a  total  of  31  configurations  are 
considered  and  up  to  six  natural  node  shapes  and  frequencies  are  given  for  each. 
This  includes  ten  configurations  of  the  “T"  substructure  and  twenty -five  configurations 
of  the  “H“  substructure,  of  which  ten  are  asymmetric. 

TO  examine  the  fully -coupled ,  three  dimensional  free  vibration  behavior  of  the 
feed  assembly,  mast,  and  solar  assembly  ("H"  substructure),  the  100 -meter 
hoop-column  antenna  of  Ref.  4  was  considered.  The  properties  used  to  model  the  "H" 
substructure  of  the  100  meter  antenna  are  shown  in  Table  5.  The  subscripts  "v"  and 
"w",  here,  signify  in-plane  and  out -of -plane  (relative  to  the  plane  of  the  antenna) 
properties,  respectively,  ta  the  case  of  the  solar  panels,  in-plane  and  out -of -plane 
designations  refer  to  the  planes  of  the  panels  themselves.  values  of  certain 
geometrical  and  material  properties  for  the  100 -meter  antenna  model  were  arrived  at 
by  scaling  properties  of  the  IS  meter  model  described  in  Ref.  6.  As  a  consequence 
of  the  lack  of  information  on  specific  solar  panel  configurations  being  considered  in 
the  design  of  the  100-meter  hoop-column  antenna,  a  solar  aasesfely  design  disriMsert 
extensively  in  the  literature7  »••*,  was  chosen  for  the  analysis.  This  design, 
referred  to  as  the  solar  electric  propulsion  (SEP)  array  is  capable  of  generating  25 
kN  of  power.  The  SBP,  as  modelled  here,  consists  of  three  main  subcomponents  as 
shown  in  Pi g.  l.  These  are  the  two  solar  arrays  themselves,  which  each  consist  of 
triangular  extendable  lattice  structure  masts  tMdi  support  flat-fold  flexible  panels, 
the  solar  boon,  which  is  foe  canister  that  the  extendable  mast  is  stored  within  prior 
to  deployment,  and  foe  base,  which  contains  such  equipment  as  foe  solar  array 
rotational  drive.  A  description  of  the  scaling  process  used,  as  well  as  additional 
explanations  for  the  choice  of  antenna  properties  together  with  all  foe  analysis  results 
are  presented  in  Ref.  10.  The  properties  displayed  in  Table  5  represent  the  "base" 
set  from  which  parametric  variations  were  performed. 

Por  three-dimensional  analysis  of  a  symmetric  substructure,  only  one  half  of 
the  mathematical  model  need  be  analysed,  much  as  in  foe  two-dimensional  case.  The 
pertinent  three-dimensional,  half-"B”  mathematical  model  is  shots)  in  Pig.  7.  The 
appropriate  boundary  conditions  on  the  center  body  are 

for  the  symmetric  casei 

all  inplane  bending  state  variables  *  0 
all  torsion  state  variables  a  0 

for  the  antisjnmmtrlc  casei 

all  out -of -plana  bending  state  variables  *  0 
all  extension  compression  state  variables  *  0 

The  effect  of  changes  in  bending  stiffness  of  foe  feed  assembly,  equal  in  in -plan*  and 
out -of -plane  directions  on  natural  frequencies  were  calculated  and  are  plotted  Pig.  • . 
This  shows  (and  mode  shapes  confirm)  that  the  feed  assembly  stiffhsss  value  of  about 
10*  M-m*  divides  a  region  below  which  foe  first  symmetric  mode  involves  primarily 
feed  assembly  motion  and  above  tfoich  it  involves  primarily  solar  panel  motion. 


The  affect  of  solar  panel  angle  and  boos  stlifneee  variations  on  antisymmetric 
mod e  frequencies  was  predicted  and  is  shown  in  Pig.  9.  It  is  clear  that  boom 

stiffness  plays  a  significant  role  in  these  trends;  at  values  of  5.10s  tt-s2  and  below, 
the  node  shapes  show  the  solar  panels  to  behave  as  rigid  bodies  in  the  first  three 
antisymmetric  modes,  so  that  their  orientation  doesn't  affect  the  corresponding  natural 
frequencies,  primarily  involving,  as  they  do,  the  flexibility  of  the  circular  cross 

sectioned  solar  boons.  At  boon  stiffhess  values  of  3xl07  N-m*  and  above,  flexible 
notions  of  the  solar  panels  and  the  mast  begin  to  be  important,  and  their  relative 
orientations  thus  vary  these  fundamental  frequencies. 

variations  in  the  lowest  four  natural  frequencies  were  examined  as  a  result  of 
the  coupling  introduced  by  a  mass  suspended  from  the  feed-mast  Junction  but  offset 
laterally  (out -of -plane)  from  the  mast  center  line  at  varying  distances.  This  makes 
the  entire  three-dimensional  "H"  substructure  asymmetric.  The  value  of  this  mass  was 
taken  as  8*  of  the  mast  mass.  Results  are  shown  in  Pig.  10.  Namely,  only  the 
fourth  mode  frequency  was  affected .  This  overall  lack  of  impact  for  substantial 

offsets  was  unexpected.  When  it  was  influential,  ie  in  the  fourth  mode,  the  mode 
shape  indicated  this  was  through  increased  axial  tension/compression  in  the  mast 
coupling  with  out-af -plane  bending. 

In  all,  /Ref.  10  provides  up  to  four  fully -coupled,  three  -dimensional ,  free 

vibration  node  shapes  and  frequencies  for  each  of  twenty-one  different  "H" 
substructure  configurations. 

n  ( m  and  Out-of-Plane  Dynamics  of  Cable -Stiffened  Polygonal  Hoops) 

This  phase  of  the  research  investigated  application  of  the  transfer  matrix  method 
to  the  determination  of  the  free  vibration  behavior  of  a  series  of  cable -stiffened  hoop 
platforms.  Hexagonal  spacecraft  subassemblies  at  this  kind  were  studied  by  Belvin11 
both  experimentally  and  using  a  NASTRAM  analysis. 

in  this  analysis  (as  in  the  complete  development  of  Phase  A),  cable  stiffeners 
were  "divided  into  half"  and  applied  at  both  ends  of  a  hoop  segment,  so  that  the 
transfer  represented  by  the  matrix  (T)  would  be  symmetric ,  so  far  as  cable  stiffeners 
were  concerned .  Note  that  the  cable  intersection  point  is  considered  fixed  to 
"ground"  in  these  studies;  in  Ref.  11,  that  point  is  free. 

Transferring  completely  "around"  this  substructure  i.e. ,  successively  across  all 
of  the  hoop  segments,  leads  to  the  requirement  for  the  structure  to  "dose"  an  itself. 
Conpatability  of  displacements  and  equivalence  of  forces  and  moments,  lead  to  a 
matrix  form  of  the  characteristic  equation  of  the  system.  For  example,  by  first 
calculating  the  matrix  IT1  -  representing  the  transfer  from  the  stats  vector  at  one 
end  of  a  hoop  segment  (station  j)  to  that  at  the  other  (station  J+l)  •  raised  to  the 
W®*  power,  where  N  is  the  number  of  hoop  segments,  the  following  relation  may  be 
written  i 

(Slj+ii  *  CT)N(*]j 

However,  since  Cslj+tr  must  equal  (slj,  this  equation  may  be  rewritten  as 

*  0 

For  Islj  to  bo  non -aero,  the  determinant  of  {tU-lTl*)  must  vanish  for  trial 
frequencies  corresponding  to  the  natural  frequencies  of  the  hoop. 


Rather  than  a  at raitht forward  Manipulation  of  tha  determinant  of 
{CU-lTl**},  however,  we  found  it  poeeible  to  determine  the  natural  frequencies 
and  mode  shapes  of  the  whole  hoop  substructure  by  considering  the  relationship 
between  the  modal  behavior  of  the  smallest  repeating  element  (i.e.,  one  hoop 
segamnt  including  half -cables  at  either  end)  and  the  remaining  elements, 
using  Thomas'  method  of  cyclic  syotttry*2.  This  is  done  by  defining  a  complex 
hoop  joint  state  vector  and  linearly  combining  two  sets  of  state  variables 
corresponding  to  modes  which  are  orthogonal  and  which  have  the  same  natural 
frequency.  That  1st 

Czlj  =  ts)j  +  iCBlj  where  i  =  t^T 

The  periodicity  of  closed  cyclically  symmetric  bodies  makes  it  possible  to  relate  the 
state  vectors  at  the  ends  of  hoop  elements  as  follows! 

IzJj  *  e^ttJj-i 

where  y  must  assums  one  of  the  values! 
y  a  2tm/N 

where  n  is  an  integer.  The  phase  angle,  y,  may  have  (M/2+1)  possible  independent 
values  for  N  even  and  (N+l)/2  values  if  M  is  odd. 

Since,  the  transfer  relations  require  that  Izlj  -  (TMz)j-x  the  cyclic 
symmetry  relation  allows  writing 

((cosy  +  isiny)CIJ  -  (Tl)tzlj  =  0 
on  CDHzlj  3  o 

and  for  this  relation  to  have  a  non -trivial  solution,  IzijMO,  the  determinant 
of  (D1  must  vanish.  The  values  of  w  which  make  (D1  =  o  correspond  to  the 
natural  frequencies  of  the  full  hoop. 


In -plane  and  out -of -plane  natural  frequencies  were  obtained  using  this  analysis 
for  hoops  with  5  through  ll  sides.  Tables  6  and  7  show  the  properties  of  the  hoop 
analyzed  and  compere  the  present  results  and  those  of  Belvin 1 1 ,  respectively. 
Figures  lia  and  b  show  the  effect  of  altering  the  number  of  segment  sides  on  the  first 
six  in-plane  and  aut-of-plane  natural  frequencies  of  the  1  meter  radius  hoop.  The 
in-plane  frame  stiffening  brought  about  by  raising  the  number  of  segments  (and  thus 
cables)  of  the  constant  diameter  hoop  is  demonstrated  t  fairly  steady  rise  in  the  first 
six  in-plane  frequencies  can  be  observed  as  the  number  of  hoop  segments  is 
increased.  Since  this  analysis  assumss  first  order  small  dlsplacemnets,  out -of -plane 
vibrations  are  influenced  less  consistently  by  cable  stiffness. 

Both  mode  shapes  and  frequencies  were  calculated  for  5,6,8,  and  11  sided  hoops 
and  are  presented  with  full  discussion  of  the  results  in  Ref.  13.  Mote  that  torque 
and  moment  balance  about  the  cable  intersection  point  was  not  always  obvious  in  these 
modee,  (see,  for  example,  the  first  in-plane  and  out-af-plane  modes  for  the 
pentagonal  hoop  in  Figs.  12a  and  b)  but  computations  confirmed  such  equilibrium  in  all 


Comparison  of  computer  "problem  state  CFO  times"  for  the  subject  method  as 
compared  to  Belvin's  NAS  TRAM  program11,  are  estimated  to  be,  for  ten  modes  and 
frequencies,  6.1  secs  and  42.3  secs,  respectively,  for  one  configuration. 
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Phase  g  ( Investigations  of  Numerical  Difficulties) 

The  full  hoop-maypole  antenna  model  described  in  Pig.  1  and  Phase  A  of  this 
report  was  programmed  for  solution  on  an  IBM  3081  -D.  Static  solutions  were  run  to 
confirm  the  internal  equilibrium  of  cable  tension  loads  and  hoop  and  mast  compression 
loads.  These  satisfactory  checks  lent  confidence  as  to  programming  accuracy. 

Attempts  to  run  natural  frequencies ,  however ,  ran  into  numerical  difficulties . 
These  were  evidenced  by  lack  of  any  reasonable  continuity  in  the  curves  of 
determinental  values  versus  trial  frequencies!  rather  such  plots  appeared  to  be  a 
scattering  of  random  points .  Initial  checks  of  the  validity  of  inverses  within  the 
program  showed  that  the  difficulties  were  not  there.  Since  the  determinental  values 
are  customarily  calculated  as  small  differences  between  large  numbers ,  double 
precision  had  been  used  at  the  outset.  The  authors  of  Ref.  2  offered  a  modification 
to  the  trial  and  error  procedure  for  finding  natural  modes  and  frequencies,  which 
ameliorates  the  "small  differences  between  large  numbers"  problems.  Essentially  it 
consists  of  assuming  a  state  vector  at  one  of  the  boundaries  (the  "near"  end)  with 
each  assumed  value  of  trial  frequency  u,  and  carrying  a  correction  column  with  an 
unknown  coefficient  for  the  state  vector.  Iterations  then  are  made  on  the  state  vector 
using  values  of  the  correction  factor  averaged  from  those  (initially)  different  values, 
which  result  from  conditions  set  by  the  known  boundary  conditions  at  the  "far"  end. 
(See  Ref.  2  pages  204  through  213.)  This  modified  scheme  made  no  improvement  in 
the  numerical  difficulties  encountered. 

In  the  studies  of  the  isolated  hoop  structures  reported,  above,  for  Phase  D,  the 
most  likely  source  of  the  numerical  difficulties  was  revealed.  The  full  antenna 
analysis  of  Phase  A  used  the  straightforward  requirements  of  hoop  "closure"  discussed 
at  the  beginning  of  the  Phase  D  report,  above.  The  isolated  hoop  analysis  began 
that  way,  as  well. 

This  straightforward  approach  requires  formation  of  the  matrix  CT1  *,  as  noted 
earlier.  l)hrigx4  and  Davies x 5  have  pointed  out,  however,  that  calculating  the 
product  of  transfer  matrices  in  a  long  chain  of  identical  components  (in  this  case  the 
hoop  segments  with  half-cables  at  each  end)  may  lead  to  large  errors,  since  all  the 
columns  of  the  matrix  product  CT1  N  tend  to  parallelism  with  the  eigenvector  associated 
with  the  dominant  eigenvalue  of  the  component  matrix  CT1  as  N  becomes  large.  Note 
that  this  is  inherent  in  the  iteration  technique  for  calculating  characteristic  values. 

A  method  suggested  by  Davies  to  avoid  the  parallelism  phenomena,  involves  the 
creation  of  a  "supermatrix"  which  contains  blocks  of  CT1  along  its  diagonal.  This 
large  matrix  is  then  used  to  obtain  CT1  11  by  Gaussian  elimination.  While  this 
approach  may  be  essentially  free  of  numerical  difficulties  for  non -closed  periodic 
structures,  problems  still  arose,  in  the  form  of  erratic  behavior  of  the  determinant 
values,  when  this  "supermatrix"  scheme  was  applied  to  the  hoop  model.  This  suggests 
that  there  are  some  further  inherent  numerical  difficulties  present  in  the  application  of 
the  transfer  matrix  method  to  closed  cyclically  symmetric  structures,  tn  this  research 
it  was  only  possible  to  avoid  them  by  taking  advantage  of  the  symmetry  of  the 
structures,  as  demonstrated  by  the  isolated  hoop  analysis  in  Phase  D. 

Thus,  the  transfer  matrix  analysis  method  developed  in  Phase  A  for  complete 
hoop-maypole  antenna  structures  would  seen  to  require  reformulation  so  as  to  allow  the 
cyclic  symmetry  of  the  suspended  hoop  to  be  accounted  for  in  the  integrated  structure, 
as  it  was  for  the  isolated  hoop. 


In  any  event,  the  user  of  a  transfer  matrix  analysis  should  be  aware  of  the 
existence  of  asymptotes  in  the  determinantal  value  vs  trial  frequency  plots  for 
mathematical  models  involving  "branches",  which  provide  multiple  paths  of  transfer 
leading  to  the  same  point.  A  simple  example  of  this  kind  is  the  "T"  substructure 
discussed  under  Phase  c,  above.  There  are  two  paths  by  which  the  transfer 
procedure  can  reach  the  single  juncture  at  the  top  of  the  "T".  Thus,  the 
determinantal  value  will  change  sign  between  the  values  of  trial  frequency  that  are  on 
either  side  of  the  natural  frequencies  of  the  cantilever  modes  of  one  half  of  the  top 
of  the  "T".  This  change  of  sign  does  not  signify  a  natural  frequency  of  the  complete 
~T~  substructure,  however,  as  shown  in  Figure  13.  While,  at  first  glance,  a 
complication  in  seeking  modes  for  the  complete  system,  knowledge  of  the  position  of 
these  asymptotes,  by  partial  analyses  of  any  system,  can  provide  guidance  which  could 
reduce  total  running  time  for  the  modes  of  a  complete,  multibranched  system. 

Conclusion 

The  research  reported  here  has  provided  insight  regarding  the  complexity  of  the 
natural  modes  and  frequencies  of  hoop-maypole  antenna  systems  by  dealing,  in  some 
detail,  with  the  free  vibratory  characteristics  of  some  of  its  substructures  and 
simplified  versions  of  the  integrated  system.  Some  typical  designs  reveal  solar  panel 
support  booms  as  being  a  critical  structural  element  for  some  fundamental  modes,  and 
in  others  the  mass  of  feed  assemblies  is  shown  to  be  a  dominant  factor.  Cable 
stiffening  effects  are  such  that  in-plane  hoop  modes  are  generally  raised  as  the 
number  of  hoop  elements  is  increased.  For  the  kind  of  symmetry  found  in  a  typical 
design,  symmetric  and  antisymmetric  modes  can  be  expected  to  be  "paired",  ie  one  of 
each  kind  found  at  nearly  the  same  natural  frequencies. 

The  transfer  matrix  analysis  has  been  found  to  be  a  viable  alternative  to  FEN 
analyses  for  the  various  substructures  of  the  hoop-maypole  type  of  satellite  antenna. 
In  fact,  considerable  computer  running  time  and  storage  requirement  reductions  relative 
to  NASTRAN  can  be  expected.  On  the  other  hand,  a  complete  analysis  still  awmits  a 
reformulation  of  the  integrated  structure  which  takes  advantage  of  the  cyclic  symmetry 
in  the  hoop  platform,  in  order  to  overcome  the  numerical  difficulties  associated  with 
analysing  closed,  repetitive  structures  with  the  transfer  matrix  method. 
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Table  1 

two  Diaeneional  Model  Geometric  Properties 
For  Baseline,  Rigid  Suspended  Hoop,  Sye—tric  Model 

Length  ( in ) 

Mass  (lb  ) 

El  (lb-ln2) 

EA  (lb) 


Peed  Mast* 1 ) 

900 

540 

16.24X109 
60 . 15x10® 


Center  Mast*2) 
2500 
1500 

16.24X109 
60 . 15X106 


Solar  Mast(3) 
900 
540 

16.24X109 
60 . 15X106 


Solar  Array  and  reed  Assembly  Mass  Taken  as  Zero. 

Cable  spring'  Rates  *  75*/ inch 

Rigid  Boop  (Spring  Suspended  Bar)t  Maas  =  800  lb 

t  Maes  Mo— nt  of  Inertia  =  4xl09  lb  -in2 

(1)  above  cable  attach— nts 

(2)  between  cable  attachaente 

(3)  below  cable  attach— nts 


1 


Table  2 


Section 


Base  Case  for  the  Upper  "T"  Substructure  (See  Pig.  l) 

Outboard  Intermediate  Inboard  Upper* 

peed  Arm  Peed  Arm  Peed  Atm  Mast 


Mass,  Kg 

68.0 

68.0 

68.0 

175.0 

Length,  m 

1.7 

1.7 

1.7 

25.0 

El,  N^npdO *® 

1.81 

1.81 

1.81 

25.40 

EA,  KxlO'7 

217.00 

217.00 

217.00 

86.00 

Table  3 

A 

Base  Case 

for  the  Lower  "T"  Substructure  (See  Pig.  1) 

Solar 

Solar  Panel 

Solar  Panel 

Lower 

Section 

Panel 

Boom 

Base 

Mast** 

Mss,  Kg 

20.0 

5.6 

305.6 

60.0 

Length,  m 

15.5 

10.9 

10.9 

12.0 

El,  N^npclO 

.001 

.001 

.001 

13.30 

EA,  HxlO*7 

.001 

3.44 

3.44 

68.80 

Table  4 

Base  Case  for  the  "H"  Substructure 

Section 

Solar 

Panel 

Solar  Panel  Solar  Panel  Pull 

Boom  Base  Mast 

Peed 

Assembly 

Mass,  Kg 

20.0 

5.6 

305.6  644.0 

204.0 

Length,  m 

15.5 

10.9 

10.9  92.0 

5.1 

El,  N^IO'6 

.OOl 

.001 

.001  25.40  1.81 

EA.HxlO'7 

.001 

3.44 

3.44  86.00  217.00 

*  from  the 

upper  cable 

attachment 

to  feed  assembly 

**  from  the 

lower  cable 

attachment 

to  solar  panels 

Table  6  -  Hoop  segment  and  cable  properties 


Mass/lengthi 

segnent  cross-sectional  A< 
El  in-planet 
El  out -of -plane : 

EA: 

GJt 


Hoop  Segments 


.8755  kg/meter 
3.2258  x  10'4  meter2 
76.96  Newton -meter2 
4925.  Newton -meter2 
22.90  Miewton 
113.3  Newton-meter2 


Cable  Length: 

Stiffness t 

a  Preload  (in -plane  case)t 
b  Preload  ( out -of -plane ) t 


Cables 


1  meter 

1.693  x  10s  Newtons 
7.5954  Newtons* 
97.224  Newtons* 


*  Values  shown  correspond  only  to  the  hexagonal  caset 

a  in-plane  cable  preload  varies  such  that  ratio  of  resulting  compressive  force 
in  each  hoop  segment  to  its  Euler  buckling  load  remains  at  .01  as  number  of 
hoop  sides  is  varied. 

b  Out -of -plane  cable  preload  varies  such  that  ratio  of  resulting  compressive 
force  in  each  hoop  segment  to  its  Euler  buckling  load  remains  at  .002  as 
number  of  hoop  sides  is  varied. 
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Table  7  -  Comparison  between  NASTRAN  and  transfer  matrix  calculation  results  for 
the  natural  frequencies  of  the  hexagonal  hoop 


i 


Mode 

TK  method 

In-olane 

NASTRAN  (Ref.  11) 

l 

92.1 

92.1 

2 

116.5 

116.6 

3 

162.0 

a 

4 

210.0 

210.1 

Out -of -plane 

itode 

* 

TM  method 

NASTRAN  (Ref. 11) 

l 

59.00 

58.41 

2 

129.00 

128 . lO 

3 

740.00 

740.30 

4 

955.00 

947.60 

5 

1369.00 

1355.00 

*  Mot  reported  in  Be Ivin’ a  results 11 ,  probably  due  to  a  difference  in  modeling 
of  the  cable  junction  at  the  "origin"  or  center  of  the  hexagon,  which  shows 


appreciable  motion  in  this  mode. 
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FIG.  I  Cross  Section  of  Hoop-Column  Antenna 
(Through  Planes  of  Feed  Assembly  and 

Solar  Panels.)  (from  Ref.  //4) 
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Figure  #  3a 
o  Bending  Mode  1 
a  Axial  Mode  1 

*  Bending  Mode  2 

•  Bending  Mode  3 
■  Bending  Mode  4 

8 

Sr 


loo  oo  rioao  moo  oo  tzso  oo  iooo.ao 
Length  a t  Central  Column  (in.) 


Figure  #  4,  18. 

a  Bending  Mode  1 
a  Axial  Mode  1 
*  Bending  Mode  2 
«  Bending  Mode  3 
■  Bending  Mode  4 


Figure  #  3b 
-  Bending  Mode  3 
»  Bending  Mode  f 
•  Bending  Mode  7 
x  Bending  Mode  8 


Figure  #  4b 
.  Bending  Mode  3 
«  Bending  Mode  3 
•  Bending  Mode  ? 
x  Bending  Mode  8 
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NATURAL  FREQUENCY  (HZ) 
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